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Narrow energy band of the electronic spectrum in some direction in low-dimensional crystals may 
lead to a negative differential conductance and N-shaped I-V curve that results in an instability of the 
uniform stationary state. Well-known stable solution for such a system is a state with electric field 
domain. We have found a uniform stable solution in the region of negative differential conductance. 
This solution describes uniform high frequency voltage oscillations. Frequency of the oscillation is 
determined by antenna properties of the system. The results are applicable also to semiconductor 
superlattices. 

PACS numbers: 73.21. Cd, 73.21.Ac, 72.80.-r 

I. INTRODUCTION 

There is a number of materials with narrow energy band of the electronic spectrum for electrons moving in certain 
directions. The narrow band may result either from low-dimensional structure of the crystal or from artificially made 
periodic superstructure like semiconductor superlattices. Typical examples of low-dimensional crystals are layered 
materials such as high-Tc superconductors, transition metal dichalcogenides like TaS2, NbSe2 and others*. Such 
materials can be treated as a stack of conducting layers of atomic thickness, separated by insulating ones*. Narrow 
electronic bands are typical also for quasi-one-dimensional conductors in directions perpendicular to conducting chains. 
For example, the linear chain conductor NbSe3 is a strongly anisotropic crystal with the anisotropy of conductivity 
Cb/cTc ~ 10 in (b-c) plane, and <7b/<T a * ~ 10 4 in (a-b) plane at low temperatures*. Such a strong anisotropy reflects 
weak inter-chain coupling and narrow electronic bands for electron motion across conducting chains. 

This feature of energy spectrum should lead to similarities of the effects observed in low-dimensional materials and 
in semiconductor superlattices which are thoroughly studied over last few decades*. Interesting effects are expected 
if electron motion approaches a regime close to the Bloch oscillations, when applied voltage is high enough and an 
electron has enough time to perform several oscillations within the same mini-band until it is scattered, i. e., when 
period of the Bloch oscillations is smaller then the momentum scattering time, t p . In such a regime DC current 
decreases with DC voltage increasing, resulting in a region of negative differential conductance (NDC). This result 
was predicted theoretically for superlattices by Esaki and Tsu?, and later was observed repeatedly. Similar conclusions 
were made for layered conductors^. However, to our knowledge, there is no direct experimental confirmation of this 
effect in low-dimensional conductors. A possible reason for this is that the crystals were not perfect enough to 
ensure the band motion of electrons in directions of low conductivity in low-dimensional materials. However, the 
strong conductivity maximum at zero bias voltage observed recently in NbSe3 by Latyshev et ali argues in favor of 
possibility of the NDC in such materials. 

An important corollary of the NDC is an instability of uniform stationary state in the NDC region which, in 
principle, may result in generation of high frequency oscillations when the DC voltage gets into the NDC region. The 
boundary of the NDC region in the static I-V curve is determined by the voltage related to the momentum scattering 
rate r" 1 , and the upper frequency limit for the NDC is equal to t" 1 as well. The value of t" 1 in many materials falls 
into the THz region£?S, so one can expect microwave generation in this frequency range. However, there is a number 
of obstacles preventing generation of high frequency oscillations. The primary one is formation of nonuniform electric 
field distribution (field domains) which impedes high frequency oscillations due to relatively low domain velocity*. 
A number of techniques was suggested to overcome this difficulty in superlattices, including application to a sample 
additional large ACS and using of pulse voltagei*, or formation of superlattices with minibands having non-sinusoidal 
dispersion law**, etc. 

In this work we study theoretically a low-dimensional conductor with current flowing in the direction of low con- 
ductivity when applied DC voltage falls into the NDC region, and find a uniform oscillatory solution in case when 
no external time-dependent influence is imposed. This solution is possible because of interaction of a sample with 
an environment which we describe by means of an effective antenna impedance. We found conditions for existence 
of stable oscillations: the antenna reactance on main oscillation frequency must have the sign corresponding to an 
inductance, in this case the sample and the effective antenna form an effective oscillatory circuit. The frequency of the 
oscillations is determined by the antenna reactance and by material properties of the sample. Oscillation amplitude 
depends on applied voltage and vanishes at the boundaries of the NDC region. 
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The paper is organized as follows: in Sec.[n]we derive general expression for the current across the layers. In Sec. lIIII 
we study in details the uniform current flow through the sample, derive equations of motion for time-dependent 
voltage, and find a solution describing oscillations. We show first the possibility of coherent voltage oscillations 
using oversimplified but physically transparent model, and then bring strict solution in the form of harmonic series. 
In Sec. IIVI we examine response of the system to a small external time-dependent perturbation, and in Sec. 
we demonstrate that solution describing coherent voltage oscillations in the sample is stable with respect to small 
fluctuations. Finally, in Sec. IVII we discuss conditions for experimental realization of the oscillatory regime found in 
preceding sections and make conclusions. 



II. MAIN EQUATIONS 

As it was mentioned in the Introduction, electron motion across the layers of perfect enough layered conductor, or 
between the conducting chains in quasi-one-dimensional conductor can be considered as a motion in narrow conducting 
band that should have specific features typical for electron motion within the minibands in semiconductor superlattices. 
In the theory of miniband transport there are two main limiting cases. If electron wave functions in adjacent layers 
are strongly overlapped then an electron has enough time to perform several tunneling processes until it is scattered. 
This quasi-classical case of miniband conduction can be well described by means of Bolzman equation which holds 
when Tp 1 <C A, where A is the miniband width. If the interlayer hopping integral is small then the sample may be 
treated as a stack of weakly coupled quantum wells, hence, tunneling processes in adjacent junctions are practically 
independent, and purely quantum case of sequential tunneling is fulfilled. An amazing result of supperlattice theory 
is that in both these limits expressions describing the current as functional of voltage are the same 4 . Below we show 
how the similar expression can be derived for a low-dimensional conductors. For definiteness we will speak about a 
layered conductor with voltage applied in the stack direction. 

For brevity we use units with e = fi = 1. To investigate non-equilibrium transport we use the Keldysh 
technique for Green functions 1 ^. Current density between layers n and n + 1 may be expressed in terms of 
the non-diagonal component of Green's function both on temporal indices and layer numbers: j n .n+i(i~) = 
ietj.fi -1 (G^„+i(r, r) - G^ +hn (r, r)) , 

(~<R (~iK \ 

^77,771 ^77,777 \ 

C A I ' 

U ^77,777 / 

where t± is the hopping integral for electron tunneling between the adjacent layers, which is supposed to be in- 
dependent on layer number and electron momentum, and G R ^ is the retarded (advanced) Green function, r is 
the dimensionless time, which is mesured in units of t p . For convenience we introduce other dimensionless units: 
frequencies and energies we mesure in terms of t -1 . 

Usually, the hopping integral, t±, is small and it is possible to construct a perturbation theory, which allows us to 
express non-diagonal on layer number Green functions in terms of diagonal ones. In the first approximation on t± it 
reads G„ !rl +i(T, r') = t± J G n . n (T, r")G n +i !rl +i(r", t')cIt", and, hence, we get well known expression for the current 
between the adjacent layers 1 ^. 




j n ,„+i(r) - it 2 x J G« n (T,r')G 7 f +1 , n+1 (T',r) + G^ n (T,r')G^ +lin+1 (r',r)dT'. 

We consider the case when electrons inside the layer are in in the equilibrium while layers are not in the equilibrium 
with each other. A difference between j n -i,n and j n , n +i leads to deviation <5/i„ of the chemical potential in layer n 
from its equilibrium value, fi. In this case in each layer the diagonal Green function with respect to layer number, 
G^n, can be related to the retarded and advanced Green functions by the expression valid for the equilibrium 
case: G^ n {ui) = (G R n (ui) — G^ n (uj)) tanh(w — 5fj, n )/2T, here we take into account temporal uniformity, so that 
G^lf'(w,w') = Gn^n\^)S(io — uj'). The transition to temporal representation we perform according to 

G^(r'y') = J e^^^^G^\J,J')J>^ dLO ' (ko ^ 

here (r) is time-dependent scalar electric potential of corresponding layer. To give proper account to scattering, we 
use the standard expression for Green functions with non-zero momentum relaxation rate G R ^n\^) — (lo — ^pzbi/2) -1 , 
here "±" is related to the retarded (advanced) Green function, respectively, £ p is the electron energy, counted from 
the Fermi energy, p is the component of electron momentum parallel to the layers. 
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We assume for simplicity that r p and t± do not depend on electron momentum. We concentrate on the case of the 
same voltage on all junctions, taking the non-uniformity into account only in Sec. [V] where we study the stability of 
the uniform solution in respect to small perturbations. So we assume that deviations, Sfj, n , of the chemical potential 
from the equilibrium value are small. 

Finally, the expression for the tunneling current between adjacent layers reads 

3 { n\r) = y[J dTxe T ' si " / ( ' 

\0 

oo n 

(Sfin+i - 8iJL n ) J drxe~ Tl cos J (V„(t) - \ 'Jt - T 2 ) i <It 2 | . (1 ) 



where s is the lattice period in the direction of the current. The second item in square brackets in Eq. \Q describes 
the diffusion current, V n (r) is the voltage between layers n + 1 and n, ctq = me 2 t 2 L T p /Ti is conductivity at zero bias 
and frequency, and s is the period of the structure. 

Consider now limiting forms and simplifications of equation JJJ for the current. We represent the voltage across 
the layer as a sum V(r) = vq + v(r), where vq is the average voltage (here we drop index n). When the second 
time-dependent term is small Eq. Q can be expanded in series with respect to i>(t), and we obtain for the current: 

j«(r) =j + - s (a d v{r) + \a' d v{rf + ^>(r) 3 ) + (c t + c' t v(r))d T v(r). (2) 

Here we introduce the DC part of the current, jo( u ) — avo/s where a = co/(l + Vq). Further, = j' a = djo/dvo = 
cr (l — t'o)/(l + Vq) 2 is the differential conductance, and c t = (3uq ~ l)/ s (l + v o) 3 is the specific tunneling capacity 
of the sample. The apostrophe denotes derivative on voltage Vq. 

As an illustration, consider, first, the simplest case when the applied voltage is time- independent, w(r) = 0. Then 
Eq. © results in well-known form£: j = jo ~ vq(1 + ^o) -1 - It follows from this expression that the differential 
conductance becomes negative when vo > 1. 

In case of small harmonic voltage, v(t) = v\e lu)T the expression for the current in the first approximation in v(t) 
reads 

From this expression it follows that at high frequencies, w > 1, the differential conductivity becomes positive even 
when the DC bias is in the NDC region. 



III. COHERENT VOLTAGE OSCILLATIONS 



In this section we derive equation of motion for the voltage and study dynamics of the system when currents in all 
junctions are identical. In addition, we neglect voltage fluctuations here, therefore, we do not consider non-equilibrium 
deviations of chemical potential in different layers and, hence, the diffusion current. 

In addition to the tunneling current Q there are other contributions to the total current. First, we must take into 
account the displacement current 

(4) 

where we use notation for the specific capacity of the sample Co = l/4ns. 

We must take into account also leakage currents, contributions due to incoherent tunneling, and other processes 
which arc not related to coherent tunneling. To allow for these effects phenomenologically we add an extra ohmic 
term to the current: 

• (ohm) /\ _ /r\ 

/ >{t)=(Ti . (5) 

s 

For the I-V curve to have the NDC region conductivity o\ must be small enough: 



a\ < a /8. 



(6) 
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To determine two independent variables, the current and the applied voltage, an additional equation relating the 
applied voltage to the current in the external to the sample circuit is required. The equation for current balance 
which determines the voltage dynamics reads 

4. ji ohm ) _|_ j( dis P l ) _|_ ji ext ) — (7) 

Here j( ext ~> = j A + j s is the current in the circuit external to the sample. We consider this circuit as that consisting 
of a power supply and environment. Coupling of the sample to the environment we describe by means of an effective 
antenna with the impedance Za- The equivalent scheme for the system under consideration is shown in Fig. ^ Since 
typical frequencies of the problem is high, we assume that AC current does not flow through power supply, so that 
antenna current ja is AC and power supply current js is DC. 




FIG. 1: Equivalent circuit of the system. Za and Zs are impedances of the antenna and of the power supply, j, jA, and 
js are the currents flowing through the sample, antenna, and power supply, respectively. Current in the external circuit 
j(ext) _ j s _j_ assume on frequencies of the problem \Za\ <C \Zs\ and on zero frequency \Za\ 3> \Zs\ so that j'a is AC 

and js is DC. 



Sample dimensions are assumed to be small with respect to radiation wavelength and to the skin-effect length. This 
allows us to neglect spatial non-uniformity. 

The resulting equation describing dynamics of the voltage oscillations has rather complicated form of the integro- 
diffcrential equation. In order to find an analytical solution we consider the limit in which the amplitude of the voltage 
oscillations is small in comparison to the DC voltage, v(t) <C Vq. This regime occurs near the boundary of the NDC 
region or in the case when the NDC region is small. 



A. Simple model 



It is instructive to examine, first, an oversimplified, but physically transparent model, based on assumption that 
in the frequency range under consideration the effective antenna can be described by the admittance with purely 
inductive reactive part: jA{v) — > -^a J v(r)dT. Here L is the antenna inductance, and factor N is number of layers in 
the sample, so that the total voltage across the sample is equal to Nv{t). The real part of the antenna admittance 
leads to redefinition of the conductivity o\ only, therefore, we drop it here. For convenience we introduce here a new 
variable related to the uniform alternating electric field in the sample as 9 = v(t)/s, and assume that the NDC region 
is small, so that contribution of the second harmonic into the friction terms is negligible (see the next subsection) . 
Eq. J7J then reads: 

e + ^e + ^(se + ^ + ^-j + a 1 v -A=o. (8) 

c s s \ 2 s J 

Here j is DC current flowing through the sample, loq = ^/ N/LSc Sl S is the sample area in the plane perpendicular 
to the current direction, and N is the number of junctions, and we have introduced the total specific capacity of the 
sample, c s = cq + Ct, and supercriticality, 5 = ad + <J\, which changes its sign from "— " to "+" across zero when the 
DC voltage leaves the instability region of the I-V curve. 

Equation for oscillations (jHJ can be treated as a mechanical equation of motion. In this equation 9 plays a role of a 
coordinate, the last term in the left-hand side of Eq. (JHJ) corresponds to a time-dependent alternating friction. In the 
regime of stable oscillations the friction averaged over the period must be equal to zero, and its momentary value can 
be considered perturbatively. So, in the first approximation we can neglect the friction term, 4- uj^Qi 1 ) = 0. Since 
the phase of the solution is arbitrary, for definiteness, we choose it in the form — E\ siwujqt. E\ is an amplitude 
of electric field oscillations on main frequency, so voltage oscillations amplitude on this frequency is E\s. To ensure 
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that the solution does not increase with time the friction term in JSJ must not contribute at the oscillation frequency 
luq. This condition determines the amplitude of the oscillations 




(9) 



Guided by mechanical analogy mentioned above, we multiply friction term in Eq. JHJ by 9\ and integrate over the 
period and demand that the work of the friction force over the oscillation period is zero. The condition of vanishing 
of the average friction determines the TV curve: 



j(«o) =J a tat-2-V. 



(10) 



This dynamic TV curve intersect the static one j s tat = (ci + <tq/(1 + Uq)) f o at three points: at the boundaries of the 
NDC region and at vq — \/3 (see fig-EJ). 




-0.6 -0.4 -0.2 0.2 0.4 0.6 



FIG. 2: Dynamic I-V curve is shown with the bold line. The dashed line describes the static I-V curve. Current is measured 
in dimensionless units, voltage 5v is counted from wo = v3 corresponding to the middle of a small NDC region. The I-V curve 
is presented for <ti = 0.11. 



B. Solution in the form of harmonic series 



In this subsection we consider a more general case and determine conditions for existence of the oscillatory solution. 
We seek the solution of Eq. JJJ for v(r) in the form of harmonic series 

oo 

v(t) = y~]( v »,n sin(nuJor) + w C: „ cos(nu! T)), (11) 

71=1 

In the most general case the antenna current reads 

N 00 

3A = — y^A a An V s,n + °An v c,n) SinnLOT + ((7An v c,n ~ ^An v s,n) COS TILOT. 
S 1 

Here a^ n = R n s/ (R^ + X^)S, a An = X n s/ (R^+X%)S, where R n and X n are resistance and reactance of the antenna 
at frequency uujq respectively. 

Solution HI 1(1 has an arbitrary phase which can be determined by initial conditions. We may fix this phase without 
any restriction assuming that v c \ = 0. The supercriticality parameter, S n — Xa An + <Ji +<r ( ji we assume to be small at 
the main oscillation frequency: \8i\ <C 1. As it will be shown later, amplitudes of the n-th harmonic are proportional 
to <5"^ 2 , hence, it is possible to neglect the terms of the order higher than the third one in the small parameter \f5\. 
So inserting Eq. 1(11(1 into the expression for the current JU, and using equation for the current balance 10, we obtain 
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a system of linear equations for v 2 s 1; v Sl 2 and v c< 2'. 



0^,2 - c' t SLUv c , 2 = 2Na J A1 - 2c s 
4cr 2 i; s ,2 - 4<S 2 w C:2 = 0. 



(12) 



Here (72 = 2c s w — Aer^. Physically reasonable oscillation frequency is determined by the second equation of the 
system (jl"2l : 



Mi 



(13) 



Here we take into account that the contribution of the second harmonic is of the second order in \/6i, so that left- 
hand sinde of this equation should be considered as higher order correction and, hence, should be neglected in main 
approximation. For oscillation amplitudes we find 



<i = -8*1 
v. % 2 = 2<5i 



<J2 2 + SJ 

V 

oS2c' t S — Cr' d 02 



Vc.2 



-2St 



V 

2ujc' t sa 2 + 52<j' d 
V 



(14) 
(15) 
(16) 



v = -3a' d ujc t scT2 + (2uj 2 c t 2 s 2 - o d 2 )8 2 + (erf + 5 2 )a d . (17) 

This expressions have a transparent physical meaning: if T> > the existence of the oscillatory regime is determined 
by the sign of the supercriticality, S±. Oscillations are possible, (i.e., the amplitude of the first harmonic is real) only 
inside of the NDC region. 

Now we show that T> is positive under conditions of our assumption that harmonics drop rapidly. If the NDC region 
is small then the voltage is close to the center of the NDC region, \vq — \/2>\ <C 1. Then a' d ex Vq — 3 is small and V 
is positive. If the NDC region is large then our approach is valid near the boundaries of the NDC region. Near the 
low- voltage boundary Vo ~ 1. In this case a d , c' t — * as v — 1 and a' d ~ 1/2, a d ~ 3/2. So if S2 and/or 72 are large 
enough, then V > 0. The high-voltage boundary of the NDC region corresponds to V3> 1. In this case c' t , a' d ~ v^ 3 , 
a'J ~ Vq A , and two first terms in l|17fl are small in comparison to the last one and, hence, T> > 0. In other words, in 
this case the influence of second harmonic on the amplitude of the first one is small, so that v s ,i — y— 861/0%, as it 
was found in the previous subsection. 

At the end of this chapter we consider the amplitudes of the higher harmonics. Acting perturbatively, we can 
express the amplitude of the third harmonic in terms of amplitudes of the first and the second harmonics considering 
the last ones as an "external force" : 

-°i3°">M + 12<J 'd( a A3 V c,2 ~ $3V S ,2) 3/2 Q , 

Vc,3 = -v„,i — —3- Ol , (18) 

24(51 + ^3 ) 

12CT d(^3^,2 + faVcf) ~ 5 3 0dVll .3/2 n a , 

V s ,3 = -V s ,l " r-g ~ 0{ . (19) 

24(51 + ^3) 

Amplitudes of the higher harmonics, n > 3, can be found in the similar perturbative way, and one can see that their 
amplitudes are proportional to S^ 2 . Therefore, the harmonic series (lllfl converges and our perturbative approach is 
valid provided <C 1. 

Thus in this section we have found the solution describing voltage oscillations in the vicinity of NDC region 
boundaries or for whole NDC region if it is small. Beyond the NDC region boundaries the oscillatory solution does 
not exist, this is reflected by the imaginary value of the oscillation amplitude in Eq. 114(1. Amplitude of the oscillations 
vanishes when the average voltage approaches NDC region boundary. Oscillation frequency is determined by geometry 
and by material properties of the sample. This solution exists provided the imaginary part of impedance of the effective 
antenna has the sign corresponding to inductance. 
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IV. MICROWAVE RESPONSE 



In this section we determine modification of the I-V curve when the system is exposed to a weak external harmonic 
radiation of frequency uj. Again, for simplicity, we treat antenna as an inductance. Let the incident radiation induces 
voltage oscillations with the amplitude Va at the antenna output. Then the equation for the oscillating component 
of the voltage (JSJ acquires the form 



2 a d° + 6 



P 3 _J- (<7i + °2)v /s 



+ /sinwT. (20) 

Lie i5 C- o <3 

Here / = VaoJq/ojNs is the amplitude of the external force. This equation describes the forced oscillations of the 
system with single degree of freedom in case of periodic perturbation. As in the previous section we consider the 
case when the amplitude of the voltage oscillations is small and we solve equation for 9(t) perturbatively. When the 
frequency of the external field, u>, is close to the eigen frequency, ujq, one can expect the effect of frequency capture. 
As in the previous section we consider the case when the amplitude of the voltage oscillations, a, is small. Then in 
the first approximation we seek the solution in the form of harmonic oscillations at the frequency of the external force 
but with unknown time-independent phase shift ip; 6{t) = asin(wr + ifi). Then we find from Eq. I|20|) 



ay/A 2 s + A 2 sin(uT + ip + ip) = f sin(wr), (21) 

where A s = uj 2 — uj 2 , A c = uj(S + a 2 uj 2 a'^/8)/c s s, and tanip = A c /A s . Eq. (|21|l determines the amplitude of oscillations 
as function of frequency, oj, and bias voltage, vq. We seek solution of Eq. lt2"T|l for the amplitude in the form 

a=^-(l + y), 

UJ 

where the first term describes the voltage oscillations with the amplitude of the natural oscillations, A\, given by 
Eq. 10 . For small enough values of / we find two solutions describing corrections to the oscillation amplitude induced 
by the external perturbation 




„( Wj ^) = ±^/(^) -<^-^>- (22) 



This expression is valid provided the external perturbation is small enough, |/| <C S 3 ^ 2 . 
Corresponding expressions for the phase read 



*-±~**M*r-t%-4r. (23) 

\UJ 2 -(J$\ 

Using the approach of the next section one can show that both solutions found above are stable. 

For uj = ujq the phases of the solutions are equal to ±tt/2. In this case the first solution describes voltage oscillations 
in phase with the driving voltage, the oscillation amplitude being increased by the external perturbation, while the 
second solution corresponds to the out of phase oscillations with the amplitude decreased by the external force. 

According to Eq. (|22|l the regime of frequency capture exists in the frequency range 

\uj 2 -uj 2 \<^. (24) 
Variation of the I-V curve due to irradiation as function of frequency and DC bias is described by expression 

A.? = d 2 1 y. 

Thus we can conclude that the external radiation induces two branches at the I-V curve that merge near the boundary 
of the regime of the frequency capture. Maximum variation of the amplitude and, hence, the largest shift of the I-V 
curve occur at uj = ujq- They are linear in driving voltage. 

At larger difference between the driving frequency and the eigen frequency, that is outside the regime of the 
frequency capture, corrections to the I-V curve induced by the external perturbation are quadratic in /. 
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V. STABILITY 



Now wc study the stability of the solution found in preceding section with respect to non-uniform perturbations of 
the voltage. Consider a sample consisting of N layers. As it was stated above, we assume non-equilibrium processes to 
be slow enough in the scale of the energy relaxation time, so that distribution function of electrons inside conducting 
layers does not change. Then the in-layer electron distribution can be described by the Fermi distribution function 
with energy shifted by a non-equilibrium correction to the chemical potential induced by difference of the currents 
through adjacent conducting layers. At first we relate the shift of the chemical potential to the voltages in adjacent 
layers. Taking into account quasi-2D character of the layers, we obtain the expression for the difference of chemical 
potentials on the adjacent layers 



<W( r ) 



■CMt)-K-i(t)). 



Ames 

This allows us to express current in terms of voltages across different junctions only: 

t! 



<?0 
S 



d,T\e 



Ams 



drie 



(V n {T)-V n (T-T 2 ))dT 2 

Tl 

cos J (V n {T) - V n (r - t 2 )) dr 2 
o 



(25) 



(26) 



where d^V n — V n +i + V n -i — 2V n is a discrete version of the second derivative on layer number. 

Since sample dimensions are assumed to be small, the current in the external circuit is determined only by the total 
voltage across the sample: 



j(ext] 



(r) 



(27) 



Using Eqs. (|2l>|) and (|2*7|) in Eq. J7J) we obtain the equation describing voltage dynamics taking into account dependence 
of the voltage on the layer number: 



<7A 



(28) 



We simplify Eq. I|28|) assuming that in the considered frequency range the antenna impedance can be modeled by 
pure inductance, 



OA 



Then each equation in Eq. 128|) acquire a form of Eq. JSJl 



k 



1 

C S S 



J dPn 



- j + (<ri + Na A ) 



(29) 



To show that solution describing coherent voltage oscillations is stable, we substitute 9 n (r) = 0^°^(r) + (5„(r) into Q2'jp. 
where Q(°'(t) is the solution found in Sec. IIIII and then linearize the resulting equation in 5 n . Temporal evolution of 
small perturbation 5 n (r) is described by equation: 



S n + (ff(r) 



Ddl) S n 



N ^ k 

k 



0. 



(30) 



Here <t(t) = \Ei<Jd cos(o;o''" + fo) + \ a d^i{^ + 2cos2(a;o'?" + ifio)) /c s plays the role of time-dependent friction coeffi- 
cient, and D = <jq / 4mc s s 2 (1 + Vq) is diffusion coefficient. 
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The general solution of linear equations l|30[) can be divided into a uniform and a non-uniform parts: 5 n = 5^ +S^ n ' , 
n = 1..N. The uniform part of the solution, S^°\ can be found as the general solution of the equation 

6 {0 \t) + ff(r)i^(r) +wg«W(r) = 0, (31) 

which follows from Eq. I|3UII when all 8 n are equal. The rest of Eq. (|30|l form N identical equations for the non-uniform 
part of the perturbation, 5^ n > 

S^(r) + [a(r)-Dd 2 n }S^(T) = 0, (32) 
which satisfy the condition which ensures that the total voltage over all layers is described by 5^ 

N 

£$W(r)=0. (33) 

n=l 

So we examine the temporal evolution of the total voltage across the sample and voltage fluctuations at different 
junctions separately, by Eqs. I|31|) and (|32fl . respectively. 

At first we study non-uniform perturbations described by Eq. I|32[l . Its solution can be found easily by means of 
the Fourier transformation on layer number: 

S k (r) = Cl (k) + c 2 (k) [ e-I^ T)+Dp)dT dT, (34) 



where k = 2sinfc/2, and ci(fc) and 02(h) are constants of integration. Constants c\(h) are irrelevant since they are not 
related to voltage fluctuations since voltage is determined by time derivative of the phase 6 n , so we can equate these 
constants to zero. Since the diffusive term promotes the stability (D > 0), we concentrate on the case of small h — > 0. 
Then the stability of the solution is determined by the damping term <j(t) which can be presented as <j(t) = <Jo + <j\(t) 
where cti(t) is a sum of harmonic terms, and the average value of the damping, <7q, is positive. Since it is positive, 
Fourier transforms of the perturbation (|34() decrease exponentially at r — * 00 for all values of h. 

Note that the results on stability with respect to non-uniform oscillations remain valid in the general case, when the 
external antenna cannot be described by pure inductance. Indeed, in case of antenna described by a linear operator of 
the total voltage across the sample, the uniform and non-uniform fluctuations can be examined separately similar to 
the way it was described above. Thus non-uniform perturbations cannot lead to instability of the oscillatory solution 
and the stability of the oscillatory regime is determined only by the stability of uniform oscillations. 

Now we discuss the uniform perturbation. The general solution of lj31|l is 

(r) = C iV (t) + C 2 v(t) [ -=^re- / a{T)dT dr. (35) 
J v z (t) 

Here C\ and C2 are arbitrary constants. As, by definition, 9^(t) — v(t), the first term in (|35|l describes the 
solution equal to the unperturbed solution that is shifted by time C\. This term is related to the different choice of 
initial conditions, and it is not related to stability. The second term decreases exponentially. This can be shown by 
direct calculations and demonstrated by means of the Ostrogradsky-Liouville formula relating the Wronskian to the 
coefficients of the linear differential equation, 



v(t) S(t) 
v(t) S(t) 



Ce 



-J a(T)dT 



where C is an arbitrary constant. With r increase perturbed solution exponentially converges to non-perturbed one 
which may have some phase shift with respect to original solution, i.e., original uniform solution is stable. 



VI. DISCUSSION 



In this section we discuss some aspects of experimental realization of the oscillatory regime. The condition for 
voltage oscillation derived in the previous section implies that the supercriticality must be negative, <5i < 0. This is 
possible when additional conductivities due to non-coherent tunneling and antenna do not compensate the NDC due 
to coherent interlayer tunneling, of. Eq. © 



a i + (J ai( uj o) < 00/8. 



(36) 
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Again we consider Eq. (jSJ. 
main frequency the "friction" 



in which we explicitly show the antenna impedance. As it was mentioned above, on 
term vanishes, and the equation defining the frequency of the voltage oscillation reads 



<j t p c s v + —9 \Z A X \ v = 0. 

Since specific tunneling capacity c s does not change dramatically with vq changes, all conclusions remain valid for the 
case of boundaries of large NDC region either. 

The oscillation frequency can be presented in the form 



TpC s S Zj a 



(37) 



Here Rq is resistance of the sample. 

Note also, that from equation for the tunneling current it follows that the oscillation frequency can not be large 
compared to the momentum scattering rate 1/t p . This becomes clear, in particular, from Eq. © demonstrating that 
the differential conductivity becomes positive at high frequencies, u> > r" 1 . 

According to Eq. I|37l) oscillation frequency depends on the input antenna impedance Za- Dependence Za(jjj) can 
be very diverse for different antenna configurations, so to illustrate possibility of the proposed oscillatory regime 
we consider the most simple case of the dipole antenna consisting of thin straight wire of length I a and radius 
r. This antenna can be described by the long line theoryii Then antenna wave impedance in the SI units reads 
p = ^-(ln 2 ^- - l)[fi] = 120(ln^- - l)[f2], where k = 2nl A /X, where A is the wave length. For the sending-end 
impedance this theory gives: 



Za 



Rap 2 



p 3 sin 2kl A 



R\ 



p 2 sin^ kl a 2 (R\ + p 2 sin^ kl A ) 



(38) 



Here, the dipole resistance Ra is a sum of the radiation resistance related to the antinode of the current, and of the 
ohmic resistance of the antenna's material, the latter is assumed to be much smaller than the former one. The radiation 
resistance is given by expression R% = P/I 2 = 30(2 (7 + In 2^;^ — Ci 2Ma) + (7 + In kl A + Ci AMa — 2 Ci 2Ma) cos 2Ha + 
(Si4A;^4 — 2 Ci2fc^) sin2fc^)[ri] (here P is radiation power, 7 is Euler constant, Ci and Si are integral cosine and 
sine, respectively). 

According to Eq. I|38|l the imaginary part of the antenna impedance is oscillatory function of frequency and on 
the length of the antenna arms I a- The dependence of real and imaginary parts of Z^ 1 on I a is presented in Fig.|3 
for several sets of antenna parameters. When the ratio Ia/X is in the range from 0.5 to 1, QZa is positive and the 
oscillatory solution is possible. 
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FIG. 3: Dependence of Z A X on ratio Ia/X for simple case of dipole antenna. Solid and dashed lines correspond to active 
resistance and reactive impedance, respectively. 

To estimate resistivity of the sample, we use conductivity data for NbSe3 adduced in Ref0 for NbSe3, <tq = 
O.lmOcm, and assume the sample dimensions to be 10pm x 10pm x 300A with the smallest length in the current 
direction. Then we obtain Rq = 3017. Typical values of ?RZa for the most interesting case of antenna with relatively 
large wave impedance which corresponds to high Q- value, are of the order of several k£l (see Fig. hence, the 
antenna impedance is larger then the sample resistance (Rq^Z^ 1 < 1). So conductivities a An and cr^ n are small, as 
it was supposed in previous sections. 
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For practical applications the higher frequencies are more interesting. We see several possible ways to increase 
the oscillation frequency. The first one is to increase Rq by increasing the ratio sN/S, either by changing a sample 
geometry or by increasing the number of layers. Implementation of this idea is limited by the condition for the NDC 
region to exist due to an increase of relative dissipation by antenna caused by radiation as a^i K s -^/ Another 
way is to maximize in order to optimize the ratio of the antenna arm length and the wavelength, or to use more 

complicated antenna, than the simple dipole one. In this case the limitations are imposed by real part of antenna 
impedance related to radiation power Jfe - practically interesting case of larger radiation leads to larger Jfe while 
">:' 2. 

In previous sections we assumed that the current density does not depend on the in-plane coordinate. In addition 
to the non-uniformity related to fluctuations of charge density considered in Sec. there is another non-uniformity 
which is related to the skin effect. So, strictly speaking, our results are applicable provided the sample width is smaller, 
than the skin-effect length, cj ^/2t:{o\ + cto/(1 + Vq))uq. This condition can be easily satisfied in the tera-Herz region. 

In conclusion, we have studied the problem of the current flow through the sample with narrow energy band 
when the average voltage corresponds to the NDC region. We have found analytically the stable uniform oscillating 
solution describing coherent voltage oscillations in the sample. The amplitude and the frequency of the oscillations 
are determined by dimensions, geometry, and material properties of the sample, and depend on the antenna properties 
of the system. The maximum value of the oscillation frequency is limited by the momentum relaxation rate. 
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